Prethermalization in quenched spinor condensates 
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Motivated by recent experiments, we consider the dynamics of spin-one spinor condensates after 
a quantum quench from the polar to ferromagnetic state from varying the quadratic Zeeman field q. 
We apply the Truncated Wigner Approximation (TWA) to the spinor system, including all spatial 
and spin degrees of freedom. For short times, we find full agreement with the linearized Bogoliubov 
analysis. For longer times, where the Bogoliubov theory fails, we find that the system reaches a 
quasi-steady prethermalized state. We compute the Bogoliubov spectrum about the ferromagnetic 
state with general q and show that the resulting finite temperature correlation functions grossly 
disagree with the full TWA results, thus indicating that the system does not thermalize even over 
very long time scales. Finally we show that the absence of thermalization over realistic time scales 
is consistent with calculations of Landau damping rates of excitations in the finite-temperature 
condensate. 



I. INTRODUCTION 

Advances in the field of ultracold atomic gases have 
spurred great interest in the non-equilibrium dynamics 
of quantum many-body systems. The ability to engi- 
neer paradigmatic model Hamiltonians, the near perfect 
isolation from the environment and the experimentally 
accessible timescales of evolution have made it possible 
to address fundamental questions about the dynamics of 
closed, interacting quantum many-body systems. 

Of particular interest in this context is the study of a 
'quench' across a quantum phase transition. Here, one or 
more parameters of the Hamiltonian are changed rapidly 
resulting in a non-equilibrium evolution of the quantum 
system towards the establishment of long range order. 
This evolution is accompanied by a spatially inhomoge- 
neous symmetry breaking and the formation of topologi- 
cal defects seeded by the quench. Accurate, time-resolved 
studies of such quenches are of fundamental importance 
to a wide range of non-equilibrium phase transitions. 

Recently, an instance of such a quench across a phase 
transition has been experimentally realized with quan- 
tum degenerate spin-1 Bose gases of 87 Rb [TJ[2] . At low 
magnetic fields, these multicomponent fluids are char- 
acterized by a contact interaction that favors a ferro- 
magnetic phase. At large external magnetic fields, the 
quadratic Zeeman energy (QZE) dominates the interac- 
tion and favors a paramagnetic ('polar') phase. As shown 
in Fig(l), these two phases are separated by a continu- 
ous phase transition. In the experiment, in situ images 
of the spin textures following a quench of the degener- 
ate gas into the ferromagnetic phase revealed the inho- 
mogeneous growth of transversely magnetized domains 
accompanied by the sporadic observation of topological 
defects that were characterized as spin vortices. 

Motivated by this experiment, here we consider the 
evolution of spin textures following a quench from the 
polar phase to the ferromagnetic phase. While the short- 



time growth of magnetization during such a quench has 
been analyzed [5HTT] , here we consider the evolution over 
much longer periods and study the manner in which the 
spin degrees of freedom thermalize following the quench. 
Applying the truncated Wigner approximation (TWA) 
to this problem (for an overview of the TWA method 
and its applications, see Refs. [12 US]), we analyze the 
long time dynamics of the spin degrees of freedom. 

We find that unless the QZE is quenched to the imme- 
diate vicinity of the critical value corresponding to the 
transition to the ferromagnetic phase, the system reaches 
a slowly evolving quasi-steady state with exponentially 
decaying correlations. This correlation length very slowly 
increases in time in a manner consistent with coarsening 
dynamics |14j . We also find and interpret previously un- 
explored physics of the dynamics of spinor condensates 
such as the emergence of longitudinal magnetization. 

II. THE HAMILTONIAN AND MEAN-FIELD 
PHASES 

Following the experimental situation, we consider the 
parameter regime characteristic of F = 1 spinor conden- 
sates of 87 Rb. Also, like in the experiment, the gases 
are confined in a quasi-2D geometry wherein the spatial 
extent d y of the condensate along one dimension is less 
than the spin healing length. This condition implies that 
the spin dynamics along this axis are effectively frozen. 
The starting point is the Hamiltonian 

71 = Wo + Hint (1) 

where the free Hamiltonian is 

H = J d W (- A. v 2 + V + qffj * (2) 

and ^ = (V'i, "00 j V'-i) 7 " is a three component spinor, V is 
the trapping potential, f x ,y.z are the spin-one matrices, 
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FIG. 1. (Color Online) Zero-temperature phase diagram of a 
spin-1 Bose condensate indicating the transition from a po- 
lar phase to a ferromagnetic phase at low quadratic Zeeman 
energy (QZE). The dashed line indicates the direction of the 
quench from an initial QZE q 3> 2|c2|nn to a final QZE q — qf. 



system in the continuum, in the absence of the trapping 
potential. In the ferromagnetic regime this quadratic 
Hamiltonian gives modes with imaginary frequencies, 
which correspond to unstable, exponentially growing in 
time modes [3]. 

To quantify the spin dynamics after the quench, it is 
natural to consider transverse and longitudinal magneti- 
zation correlation functions 



G x (r,t) = F + (r)F_(0) :) 
G z (r,t) = i<: F z (r)F z (0) :) 



(5) 
(6) 



and q is the quadratic Zeeman shift. The interaction 
Hamiltonian is 



n 



int 



l 



-c n' + \c 2 F' 



(3) 



Here we have n = F = ^f\E', and the param- 

eters Co and c 2 can be expressed in terms of the s- 
wave scattering lengths as cq — (ao + 2a2) and 

c 2 = t™3m (° 2 — a o)- There is also a linear Zeeman shift 
in the system proportional to F z but this can be dropped 
since F z commutes with the Hamiltonian and thus is 
conserved. Spinor gases of 87 Rb are characterized by 
the scattering lengths (20(02) = 101.8(100.4) Bohr radii 
respectively. Because a 2 < ciq, the spin dependent in- 
teraction in these spinor gases favors the ferromagnetic 
state. Since for these parameters Cq ^> | C2 1 , fluctuations 
in the total density are suppressed. Further, due to the 
weak spin interaction strength, the interparticle separa- 
tion at typical densities is far smaller than the spin heal- 
ing length, thereby suppressing the role of quantum de- 
piction. 



III. SHORT-TIME THEORY 

As illustrated in Fig. [IJ for q > q = 2|c2|n , the 
ground state is the polar state while decreasing q be- 
low q 0l the ground state acquires a magnetic moment 
and eventually reaches the fully polarized ferromagnetic 
state at q = 0. Initially, we take q = q L to be the largest 
energy parameter in the Hamiltonian and quench to a fi- 
nal state q = qj where < qf < qg. For short times, the 
dynamics is expected to be well-described by expanding 
the Hamiltonian to quadratic order around the initial po- 
lar state [SHS], ^ = y / «o(0, 1,0) T . Under this expansion 
the Hamiltonian becomes 

k 

+c 2 n o (-0 lik V'Li _ k + i>-l,-ki>l,\S) ( 4 ) 

where we have dropped the term describing the stiff den- 
sity fluctuations. For simplicity, we have considered the 



[3J and the concomitant gain functions, Gj_(t) and G z (t), 
which are the above evaluated for r = 0. For qot/h 3> 1 
and qf = 0, the above correlations within the linearized 
theory give 



G±(t) 



G z (t) 



1 



1 



aQ0t/h 



^fSirq t/h n g 
1 1 
Mn^ti 2 (n e s ) 2 ' 



,2dat/H 



(7) 



(8) 



where £ s 



H 2 /qom is the spin coherence length. In- 



terestingly, the longitudinal gain grows with twice the 
exponent of the transverse gain. However, the longi- 
tudinal magnetization is suppressed by a large factor 
l/( n o£s) ~ 10~ 4 for the experimental parameters of 
Ref. pp. Thus, G z (t) remains much smaller than G±(t) 
up to the times where the latter saturates G±(t) ~ 1 and 
the linearized theory does not work (see Fig. (2). 



IV. TRUNCATED WIGNER SIMULATIONS 

As seen from the exponential growth of the gain func- 
tions, the above theory clearly fails once the transverse 
magnetization is of order unity. To gain a more com- 
plete understanding we turn to truncated Wigner sim- 
ulations. This technique involves propagating the full 
spinor Gross-Pitaevskii Equations (GPEs) of motion ob- 
tained by taking the classical limit of Eq. |l]) seeded with 
random initial conditions for canonical degrees of free- 
dom distributed according to the Wigner transform of 
the initial state. For our case the initial state is taken 
to be the vacuum of Bogoliubov quasiparticles for a the- 
ory expanded about the polar state. For the limiting 
case of = 00 this has the particularly simple form: 
W(tt) « <5(|^ 0i0 | 2 -no)n k 8exp [-2£ Q |^ Q , k | 2 ] , where 
the factor of 8 ensures the correct normalization of the 
Wigner function [13]. 

To propagate the wave functions governed by the 
spinor GPE, we first discretize the system effectively 
by using a lattice with the constant I chosen such that 
I <C £ s and n £ 2 ^> 1 and then use a split-operator 
method that is accurate up to cubic order in the time 
step size [T5J IHj • The TWA is an approximate method 
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FIG. 2. (Color Online) The growth of transverse (red) and 
longitudinal (blue) magnetization following a quench to qf = 
0. Inset: The characteristic time constants for the growth 
of transverse and longitudinal magnetization densities versus 
the 2D interaction strength \c2n\. The dashed lines indicate 
the predictions of the linearized Bogoliubov treatment for the 
respective parameters. 



resulting in the leading order in the expansion of classical 
dynamics in quantum fluctuations |13) . It is also known 
to be asymptotically accurate at short times and exact for 
quadratic theories. In our case we expect this method to 
be quantitatively accurate at all times. Indeed the small 
parameter of the expansion is l/(no^) <C 1 so that quan- 
tum corrections to TWA are suppressed by this factor 
and are potentially only important at very long times. 
At longer times when nonlinearities become relevant, the 
momentum modes become highly occupied and the dy- 
namics remains classical. Putting the system on a lattice 
ensures that there is no potential spurious effect from 
vacuum occupation of high-energy modes, which some- 
times impedes the validity of TWA at long times [T2] , 
We checked that the results of our simulations are insen- 
sitive to the choice of the cutoff as long as it is shorter 
than the spin coherence length. 

Shown in Fig. [2] are the TWA results for the transverse 
and longitudinal gain functions in the absence of a trap 
compared with the analytic results based upon the above 
linearized theory. At short evolution times following the 
quench, the simulations predict the exponential growth 
of both the transverse and longitudinal magnetization 
densities in a manner that is in excellent agreement with 
the linearized theory. During this period, the spin tex- 
tures are characterized by ferromagnetic domains that 
are predominantly oriented in the transverse plane. 



V. LONG-TIME BEHAVIOR 

We now move on to discuss the long-time dynamics 
of the quenched spinor condensate. As is seen in the 
gain functions in Fig. [2j the magnetization of the system 
eventually saturates and reaches a steady state. This 



motivates one to consider the prospect of thermalization 
in the system. More specifically, under the assumption 
of thermalization, the system will evolve to the ground 
state at q — qf, with the excess energy accounted for by 
a superposition of elementary excitations about this con- 
figuration. We therefore will consider the theory about 
the state at q = qf. 

We define the heating (excess energy) of the system as 

Q = (U)i - (H) ES (9) 

where the above expectation values are evaluated for the 
initial state immediately after the quench and the ground 
state for q — qf. To extract the equilibrium temperature 
this heating can be compared with the thermal energy of 
the system. It is convenient to use the following param- 
eterization of the spinor 

= \fne la (sin 77 cos 4>e 1 ^ , cos rje n , sin 77 sin <j)e~ % ^ T . 

(10) 

With this, one can see that for < qf/qo < 1 the classical 
energy is minimized for (j) = 7r/4, 7 = 0, and cos(2ry) = 
qf/qo- One then obtains that the above heating is 

Q = \Nq (1 - q f /q a ) 2 (11) 

where the details of this calculation are in Appendix [A] 
Note that the above is based on classical energy dif- 
ferences. In addition to this classical contribution to 
the heating there is also a quantum correction coming 
from the zero point fluctuations of the condensate. How- 
ever this contribution is suppressed by a large factor 
U noQ ~ 100 and thus is not important. 

The thermal energy of the final state is found from 

Q = U{T) = TV (e- n/T n^j /Z (12) 

where Z = Tr (e~^/ T ) is the partition function and the 
ground state energy of the system is set to be zero. To 
evaluate the above we develop a second Bogoliubov the- 
ory by expanding the Hamiltonian to quadratic order 
about the minimum at q = qf. Note that within the cho- 
sen parametrization one automatically avoids subtleties 
related to the lack of a long range order in two dimensions 
at finite temperatures. The quadratic theory is then self- 
consistently checked by verifying that the resulting ther- 
mal depletion is small for the computed temperature. 

The full Bogoliubov analysis of the spectrum linearized 
around the ferromagnetic minimum is rather cumber- 
some because of the phonon-magnon coupling |17j (which 
vanishes for the special points qf — and qf = qo). How- 
ever, 87 Rb has a natural separation of energy scales since 
Co 3> I C2 1 which inhibits density fluctuations. This signif- 
icantly simplifies the Bogoliubov analysis and effectively 
results in the decoupling of the stiff density fluctuation 
and the spin modes. In this limit one finds the simpli- 
fied spectrum consisting of two modes with the dispersion 
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(see Appendix [A| 



(-L) 



£k(ek + g/) 



(13) 



M) - 



\ 



(£k + go) £k 




where £k is the free particle dispersion. Note that these 
reduce to the known dispersions in the limits of qj = 
and g/ = go [IB HFJ- Furthermore, the large Co limit 
of the full dispersion [17] agrees with Eq. (13). These 
expressions are used to compute the thermal energy of 
the system. It is straightforward to verify that unless 
qf is very close to qo, 1 — g//go ^> l/\/ n o£s > the dom- 
inant contribution to the thermal energy U (T) comes 

from the quadratic dispersion of these modes: w k ~ 

uj£ ps £k- This can be justified a posteriori by com- 
paring the temperature with q . With these conditions, 
U(T) fa 7V7rT 2 /6(n ^)g . Finally, equating this to the 
heating, one finds 



T = 



2tt 



«o^(go -qf). 



(14) 



To test the hypothesis of thermalization, we use the de- 
veloped quadratic theory to evaluate finite temperature 
correlation functions (evaluated at temperature Eq. (14)) 
and compare with long-time TWA results. The most nat- 
ural correlation function to consider is Gj_ defined in Eq. 
([5]). The gapless mode, Eq. ( 13 ), which leads to algebraic 



decay of the correlation function in two dimensions has 
the dominant contribution. Within the quadratic Bogoli- 
ubov theory, the correlation function is found to be (see 
Appendix [b]) 



Gj_(r)= l-(g//gof (r/&) 



(15) 



with exponent a = -\/6/7r 3 • \/l/no£f • Because of the 
dependence of the exponent a on the small parameter 
l/ n o£s: we see that, assuming equilibrium, G±(r) will 
not decay by any appreciable amount over the relevant 
length scales of the condensate, which are typically of the 
order of tens or hundreds of the spin coherence length £ s . 
For the same reason in equilibrium the system should not 
have any vortices. 

We point out that the resulting thermal state of the 
condensate with temperature T given by Eq. (14) has 



some interesting properties. On the one hand except for a 
small vicinity of g/ ss go it is a high-temperature classical 
state characterized by a temperature much higher than 
the chemical potential. This implies that quantum de- 
pletion in this state is negligible. On the other hand, the 
smallness of the exponent a shows that this state is well 
below the Kosterlitz-Thouless (KT) transition tempera- 
ture due to the unbinding of vortices in the spin degrees 
of freedom. More specifically, the KT transition temper- 
ature scales as Tkt ~ qo n o£,s which can be compared 
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FIG. 3. (Color Online) The transverse magnetization den- 
sities for short (a) and long (b) times for q/ = over a 
70 x 70/^m region. The amplitude and orientation of the trans- 
verse magnetization are indicated by the brightness and hue 
(color wheel shown), (c) The spin correlation function for 
these same times. While the correlation function at short 
times is in excellent agreement with Bogoliubov theory, the 
decay of this function at long evolution times is much more 
rapid than predicted by Eq. (|15[) (dashed line). 



with Eq. (14). In this respect, this is a low-temperature 
state. 

We now compare these predictions to the results of 
the TWA simulations for long times. Shown in Fig. [3] 
are magnetizations for short and long times as well as 
correlation functions after the quench. Displayed are 
results for zero final quadratic Zeeman field qf = 0, 
though qualitatively very similar results occur for finite 
values of gj (for a full animation of the dynamics see 
EPAPS material [2D]). As shown in Fig.^c), the corre- 
lation function Gj_ before saturation for short times has 
the functional dependence of a Bessel function Jo(r/£ s ) 
as previously predicted [3;. For long times, the cor- 
relation function reaches a steady state and decays by 
over a factor of five from its r = value for a system 
size on the order of 50 spin coherence lengths. Such 
decay is qualitatively incompatible with the theory de- 
veloped by assuming thermalization. We thus conclude 
that these spinor condensate systems do not thermalize 
at appreciable time scales but rather reach a quasi-steady 
(prethermalized) regime which evolves anomalously slow 
in time. Such prethermalized phases were suggested ear- 
lier in the contexts of weakly interacting fermions [2"Tll22| . 
Bose supcrfluids 23, 24 , and dynamics of the early Uni- 
verse [25]. What is perhaps unexpected is that in this 
system anomalously slow relaxation occurs at very high 
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levels of the heating in the system. 

To support these results, it is useful to compare with 
the Landau damping rates of long-wavelength Bogoli- 
ubov modes as a result of scattering from thermal modes 
Such an analysis is carried out in Appendix |C| 
It is found that typical lifetimes of these modes is on 
the order of 100 seconds. This rate is considerably lower 
than those of scalar condensates due to the weak spin- 
dependent interaction and small thermal depletion. The 
lifetime estimates provide lower bounds for the thermal- 
ization time and are consistent with the numerical re- 
sults. However, such long time scales are outside of our 
numerical as well as experimental access. 

In conclusion, we have examined the long time evo- 
lution of quantum degenerate spinor gases following a 
quench to a ferromagnetic phase. Assuming a thermal- 
ized final state, we find that the finite temperature spin 
correlations are characterized by an algebraic decay over 
length scales much larger than the relevant length scales 
of the condensate. In distinct contrast, numerical sim- 
ulations based on the truncated Wigner approximation 
indicate a rapidly decaying correlation function even at 
extremely long evolution times. This inconsistency leads 
one to the conclusion that the quenched spinor conden- 
sates will not thermalize over experimentally relevant 
timescales. These results are consistent with the Landau 
damping rates. The role of dipolar interactions [21 [30] in 
the long time evolution has yet to be examined. 
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Appendix A: Bogoliubov analysis in the regime 

< q f < 2\c 2 \n 

The easiest way to perform the Bogoliubov analysis 
of spinor condensates is to use a parametrization of the 
spinors though the density-angle variables as in Eq. ( 10 1. 



In this case one avoids issues related to absence of the 
true condensation in one-dimension at finite tempera- 
tures. Also the normal modes have a very transparent 
physical meaning. The analysis of this Appendix very 
closely mimics an analysis for the description of the low- 
energy excitations of bosons in an optical lattice close to 
the superfluid-insulator transition in the effective spin- 



one representation [3TJ [35] ■ In this parametrization, the 
magnetization F z = ^*/ z $ and F z 2 = 'I'*/ 2 * are given 
by 

F z = |*i | 2 - |*_i| 2 = nsin 2 ^) cos(20), (Al) 
F z , 2 = 1*1 1 2 + |*-i| 2 = nsin 2 (r?). (A2) 

Similarly the square of the transverse magnetization 
|Fi I 2 = F 2 + F 2 reads: 



|Fj_| 2 = — sin 2 (2ry)[l + sm(20) cos(27)] 



(A3) 



With these expressions, one finds for the interaction en- 
ergy density (here we include the quadratic Zeeman term 
in the interaction energy) 



C JLn 2 +°*(F 2 
2 2 y z 



\F, I 2 ) 



qF z , 



= \ 2 + c 2 n 2 sin 2 (77) cos 2 (77) [1 + sin(20) cos(27)] 



qn sin 2 (77) 



sin 4 (77) cos 2 (20). 



(A4) 



It is straightforward to check that for C2 < and < 
q < 2\c 2 \n = q the energy is minimized when = n/4, 
7 = 0, and 77 — fj where we define 77 according to q — 
qo cos (277). There is an equivalent minimum obtained by 
gauge transformation where — > <f> + ir/2 and 7 — > 7 + 7T. 
The minimum of the interaction energy density is then 



E(q) « -g « sin 4 (77) 



4 



1 



q_ 
go 



(A5) 



Note that in this expression we ignored the zero point 
energy associated with the depletion, which is suppressed 
by a large factor 1 / y/n~oQ. 

Within the Bogoliubov approximation we need to ex- 
pand the energy around the minimum. For this purpose 
we define small deviations of the angles and 77 from 
the optimal values — 7r /4 + <50, n = fj + 5n and do a 
second order expansion in S(f>, Sn, 7, and £. Since we are 
interested in the limit Co S> c 2 the density fluctuations 
are suppressed and we can set Sn = and fix the density 
at its equilibrium value n = rig. Then the expression for 



the interaction energy density ( A4 ) reduces to 



SHiat ~ 9o^o sin 2 (77) 

x [cos(2t7)((50) 2 + cos 2 (t7)7 2 + 4 cos 2 (77) ((J??) 2 ] . (A6) 

Under the same approximation the kinetic energy density 
term reads 



%kii 



2m 



(VSi]) 2 + sm 2 (fj)(V(/>) 2 + sin 2 (77)(V£) 2 



-cos 2 (t?)(V7) 2 + (Va) 2 + 2VaV7COS 2 (7?) 



(A7) 



And finally to determine the canonically conjugate vari- 
ables we need to write the Berry phase term, 7\E^* in the 
linearized approximation: 

2no sin 2 (77)50 ^ + ?7-o sin(2?7)£777 — 5n[o!-|- cos 2 (77)7]. (A8) 



G 



In this form it is clear that the phase conjugate to the 
density fluctuations is a — a+cos 2 (fj)j. Then the kinetic 
energy density can be rewritten as 



2m 



n 



(Vn) 2 + sin 2 (?7)(V0) 2 + sin 2 (77)(V£) 2 

(A9) 



+ (V<5) 2 + cos 2 (r?) sin 2 (?7)(V7) 2 

One can further redefine variables 

5<t> -> 0/(2 v / n o sin(^)), £^ £ / y/n^ sin(fj) , 



7 ^ 7/v / ^o sin ( 2 "); V^v/V^o 



Then, the spin part of the Lagrangian density (Berry 
phase minus Hamiltonian) becomes 



£ w <j> £ + ffi - g s 



2m 



(Vry) 2 + -(V0) 2 + (V£) 2 + cos 2 (fj) (V7) 2 



(A10) 



This Lagrangian density immediately gives two sets of 
normal modes corresponding to oscillations in the cj) — £ 
variables, which represent uniform rotations around z- 
axis, and in the rj — 7 variables, which represent oscilla- 
tions in the magnitude of the transverse magnetization. 
The corresponding frequencies in the momentum space 
are 



A) 



= V £ k(£k + <7ocos(277)), 



(AH) 

dj? - \J{9o + e k )(e k + q sm 2 (2fj)) (A12) 

where e k = is the free particle dispersion. We note 
that these expressions follow from taking cq — > 00 limit 
of the general dispersion relations obtained in Ref. [T7] , 
There are, however, advantages in using the density-angle 
representation since it does not rely on the assumption 
of the existence of the true long range order. 



Appendix B: Finite temperature correlation 
functions within the Bogoliubov theory 

The Bogoliubov theory described in the previous sec- 
tion makes it possible to compute correlation functions 
in two-dimensional systems. In the following we will con- 
sider the transverse magnetization correlation function 



G x (r,t) 



1 



F + (r)F_(0) 



(Bl) 



According to the Mermin- Wagner theorem, a gapless 
mode will lead to a power law behavior of correlation 
functions in two-dimensions. The gapless mode for our 
case corresponds to the 4> — £ variables. It is straightfor- 
ward to see that these will dominate the correlation func- 
tions. Taking into account these two fluctuating fields 
one finds 

< F+(r)F_(0) >= sin 2 (277)e- <A ^ 2> / 2 e- <A « 2> / 2 . 

(B3) 

In this equation, (A£ 2 ) = ((£(r) — £(0)) 2 ) with a similar 
expression for (A(j> 2 ). 

Using the analysis from above, one finds 



<A£ 2 ) = 
and 

(A^ 2 ) = 



1 



(-L) 



N sin 2 (fj) £7? 



E^T) (/(4 ±) ) + l/2) (l-Jb(fcr)) 



/(4 ±} ) + l/2) (l-Jo(Ar)). 

(B4) 



Nsm 2 (fj)^J^ 



where / is the Bose distribution function. It can be seen 
that the dominant contribution comes from (A£ 2 ) due to 
the long-wavelength divergence of the sum over k which 
is cut off at ~ 1/r. The above expression is therefore 
well-approximated by 



<A£ 2 ) 



T 

irsm 2 (fj)q n £,'; 



log(r/6). 



(B6) 



Using the expression for the temperature derived in the 
manuscript, one arrives at Eq. (15 1. 



Appendix C: Landau Damping Rate 

We take a spinor condensate with a single Bogoliubov 
excitation of wave vector k and evaluate its lifetime r^. 
due to scattering off of short-wavelength thermal modes. 
For simplicity we will consider the gapless spin mode w k 
as given in Eq. (13). Such a rate is given by the well- 
known Landau damping formula [26 which has been gen- 
eralized to Bose-Einstein Condensates in 



f = iEi^i 2 (/(4 ±) )-/(4 ± + l))^(4 ±) +^ , -4U' 

k k' 

(CI) 

where / is the Bose-Einstein distribution function. In 
this equation, the matrix element Mkk' is given by |27l - 
[29] 



Written in terms of the variables introduced in the pre- 
vious section we have 



F+ = -^ne ia sin(2n)e~^ (cos^e 47 + sin(0) e - 17 ) 
\J2 



(B2) 



,(-L) 



M kk , = J* = J*L(I* + ^ . J (C2) 

where s is the sound speed of the mode w k , . With 
this expression, the two-dimensional k'-summation can 
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be performed. In the limit T ^> qf, the result is 
1 _ a s k T q f ^ _ fcg s T q f . 



pression for the temperature given in Eq. ( 14) and taking 
qj = qo/2 we find that 



Tk 



ft d y /£ s q 



1.13 



(C3) 



where a s = (a2 — ao)/3 is the spin-dependent scattering 
length. Using the experimental parameters, and the ex- 



r fc » 90s 



for a mode having wave vector q = l/(2£ s 



(C4) 
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